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Abstract 

In this paper an asymptotic homogenization method for the analysis of composite ma¬ 
terials with periodic microstructure in presence of thermodiffusion is described. Appropriate 
down-scaling relations correlating the microscopic fields to the macroscopic displacements, tem¬ 
perature and chemical potential are introduced. The effects of the material inhomogeneities 
are described by perturbation functions derived from the solution of recursive cell problems. 
Exact expressions for the overall elastic and thermodiffusive constants of the equivalent first 
order thermodiffusive continuum are derived. The proposed approach is applied to the case of 
a two-dimensional bi-phase orthotropic layered material, where the effective elastic and ther¬ 
modiffusive properties can be determined analytically. Considering this illustrative example 
and assuming periodic body forces, heat and mass sources acting on the medium, the solution 
performed by the first order homogenization approach is compared with the numerical results 
obtained by the heterogeneous model. 

Keywords: Periodic microstructure, Asymptotic homogenization, Thermodiffusion, Overall 
material properties. 


1 Introduction 


Composite materials are extensively used in industrial practice. Indeed, many advanced engineer¬ 
ing applications, such as aerospace, aircraft, green building, biomedical, energetics and electronics 
require the design and the use of heterogeneous multiphase materials. Due to the microstructural 
effects as well as the interaction between their constituents, these materials may present several 
favorable physical properties, as for example high stiffness, improved strength and toughness, en¬ 
hanced thermal conductivity, mass diffusivity or electrical permittivity. 

Recently, multiphase composite materials have been largely used in the design and fabrication 
of batt ery devic es, in particular of lithium-ion batteries and solid oxide fuel cells (jNakaio et al 


of battery devices, m particular o t fitm um-ion batteries and solid oxide tuef ceffs I IN a ka 10 et af.l . 
2012t [Dev et all 2014; Elli s et al. . 2012 1. Since high operational temperatures can be reached and 
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intense particle fluxes are needed for maintaining the electrical current, the components of such 
battery devices are subject to severe thermomechanical stresses as well as stresses induced by the 
particle diffusion, which can cause damage and crack formation, compromising the performa nce o f 
the devices in terms of power generation and energy conversion efficiency ( Atkinson and Sunl . 20071 : 


Delette et all 20131 ). Modelling the mechanical and thermodiffusive properties of the components 


of such battery devices represent a crucial issue in order to predict these phenomena and then to 
ensure the successful manufacture and the reliability of the systems. 

The macroscopic behavior of thermodiffusive composite materials used for realizing lithium-ion 
batteries and solid oxide fuel cells is influenced by multiphysics phenomena occurring at scale-lengths 
characteristic of the m icroscopic constituents, which is small compa r ed to t he macroscopic dimensio n 
(i.e. structural size) I Richardson et al. . 20121 : Bove and Ubertini . 20081 : Haiimolana et all l201ll) . 


Consequently, multiscale techniques represent an appropriate and powerful tool for modelling the 
effects of the microstructures on the macroscopic mechanical and thermodiffusive properties of these 
materials. In particular, for composites with periodic microstructures, homogenization techniques 
represent an useful and advantageous method for providing a rigorous and synthetic description of 
the effects of the microscopic phases on the overall properties of the materials. The application 
of these approaches makes possible to avoid the challenging numerical computations required by 
computational modelling of heterogeneous media. 

Several homogenization techniques have been proposed for studying overall static and dynamic 
elastic propertie s of composite materials with periodic microstructu r es, su c h as the asymptotic 


(19891); Allaire (1992); Boutin and Auriault 1993b; Meguid and Kalamkarov 

(19941); Boutin (199fil): 

Andrianov et al. (2008j): Tran et al. (2012)1, the variational-asvmptotic methods (see for example 

Smvshlvaev and Cherednichenko 2000J); Peerlings and Fleck 20041); Smvshlvaev (20091): Bacigalupo 

(20141); Bacigalupo and Gambarotta (20141)) and the computational 

approaches (see for example 

Forest and Sab (19981): [Forest (2002); Kouznetsova et al. (2002, 2004 

); Kaczmarczvk et al. (20081); 

Forest and Trinlr (20111); Bacigalupo and Gambarotta (2010l 2011, 

2013); 

De Beilis and Addessi 


ered heterogeneous material at the micro-scale, described by a standard Cauchy continuum, an 
equivalent homogenous medium at the macro-scale. The behavior of the equivalent macroscopic 
material can be described by means of a first order continuum or alternatively a non-local medium. 
Multiscale asymptotic and computational homogenization procedures have been also proposed for 
the a nalysis of heterogeneous media in presence of multiphysic s phen omena, such as thermomechan¬ 
ical ( Kanoute et all 2008: Zhang et all . 120071 : lAboudi et all 2001 ) and thermo-magneto-electro- 


elastic ( Sixto-Camacho et all 20131) deformations. Recently, these methods have been applied for 


studying the influence of the microstructural effec ts o n t he macro s copic al mechanical behavior and 
operative performances of lithium-ion batteries ( Salvadori et all l2014h and solid oxide fuel cells 


(Bacigalupo et al.,|201J). The overall properties of periodic multilayered structures characterizing 


such energy devices can be efficiently described by means of homogenization methods developed 
for periodic composite materials. Nevertheless, to the author’s knowledge, a rigorous asymptotic 
procedure accounting for the effects of the microstructures on both macroscopic elastic and ther¬ 
modiffusive properties of composite materials as well as on the coupling between these properties 
is still unknown in literature. 

In this paper, an original asymptotic homogenization method for modelling the static elastic, 
thermal and diffusive properti es of periodic thermodiffusive composite materials is proposed. The 
rigorous approach developed in I Bakhvalov and Panasenko ( 1984 )): Smvshlvaev and Cherednichenko 
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(2000); B acigalupol (2014) and Bacigalupo and Gambarotta ( 2014 ) is extended in order to account 
for the effects of the microstructures on the macroscopic temperature and chemical potential of 
the materials and on the stresses induced by these fields. The displacements, temperature and 
chemical potential at the micro- and macro-scale are related through an asymptotic expansion of 
the microscopic fields in terms of characteristic size e of the microstructure. This expansion depends 
both on the macroscopic strains, temperature and chemical potential gradients and on unknown 
perturbation functions accounting for the effects of the heterogeneities. Perturbation functions 
representing the effects of the material microstructures on the displacement, temperature, chemical 
potential and on the coupling effects between these fields are introduced. These perturbation 
functions, depending only on the properties of the microstructure, are obtained through the solution 
of non-homogeneous problems on the cell with periodic boundary condition s. _ 

Sim ilarly to the procedure proposed in ISmvshlvaev and Cherednichenko ( 2000h and Bacigalupol 
( 2014h . averaged field equations of infinite order are obtained, and their formal solution is performed 
by representing the macroscopic displacements, temperature and chemical potential in terms of 
power series. Field equation for the homogenized first order thermodiffusive continuum are derived, 
and exact expressions for the overall elastic and thermodiffusive constants of this equivalent medium 
are obtained. The proposed formulation is applied to the case of a two-dimensional bi-phase or¬ 
thotropic layered material. The effective elastic and thermodiffusive constants corresponding to 
this example are determined analytically using the general expressions derived by the homogeniza¬ 
tion procedure. The solution performed by the proposed approach is compared with the numerical 
results obtained by the heterogeneous model assuming periodic body forces, heat and mass sources 
acting on the considered bi-phase layered composite. 

The article is organized as follows: in Section 2 the geometry of the considered thermodiffusive 
composite material with periodic microstructure is illustrated, and the corresponding constitutive 
relations and balance equations are introduced. The developed multiscale asymptotic homogeniza¬ 
tion technique is described in Section 3, based on down-scaling relations correlating the microscopic 
fields to the macroscopic displacements, temperature and chemical potential. The unknown per¬ 
turbation functions describing the effects of the material heterogeneities are defined as solutions 
of the corresponding non-homogeneous cell problems. In the same Section, averaged field equa¬ 
tions of infinite order are obtained, and a solution scheme based on asymptotic expansion of the 
macroscopic displacements, temperature and chemical potential field is reported. Field equations 
and explicit expressions for the overall elastic and thermodiffusive constants of the equivalent first 
order homogeneous continuum are derived in Section 4. As just anticipated, the proposed approach 
is applied for studying overall properties of two-dimensional bi-phase orthotropic layered materials 
in Section 5. Finally, a critical discussion about the obtained results is reported together with 
conclusions and future perspectives in Section 6 . 


2 Governing equations of periodic multiphase materials in 
presence of thermodiffusion 

Let us consider an heterogeneous composite material having periodic micro-structure and subject 
to stresses induced by temperature changes, mass diffusion and body forces. The two-dimensional 
geometry shown in Fig. [T] is assumed for the system. Considering small strains approximation, 
the constituent elements of the medium are modelled as a linear thermodiffusive elastic Cauchy 
continua. The material point is identified by position vector x = x\e\ + 2:262 referred to a system of 


3 






















v, =^e. 


Figure 1: (a) Heterogeneous material - Periodic domain L ; (b) Periodic cell A and periodicity 
vectors. 


coordinates with origin at point O and orthogonal base {ei, e 2 }. The periodic cell A = [0, e] x [0, 5e\ 
with characteristic size e is illustrated in Fig. [lj). The entire periodic medium can be obtained 
spanning the cell A by the two orthogonal vectors tq = d\ei = eei, v 2 = d 2 e 2 = dee^- 

According to the periodicity of the material, A is the elementary cell period of the elasticity 
tensor 

+ = C (m ’ £) (a;), *=1,2, Vx£A, (1) 

where the superscript m stands for microscopic field. Similarly, the heat conduction tensor K^ m,e ^ ( x) 
and the thermal dilatation tensor a ^ m,e ) (x) are dehned as follows 

K {m ’ £ \x + v x ) = K {m ’ £ \x), a^ m ' £ \x + v i ) = a^ m ’ £ \x), * = 1,2, Va: e A, (2) 

and then the mass diffusion tensor and diffusive expansion tensor /3^ m,e \x) become 

D {m ’ e \x + v,) =D {m ’ £ \x), P {m ’ £ \x + v i )=f3( m ’ e \x), * = 1,2, Va; e A. (3) 


The tensors ©, © and © are commonly referred to as A—periodic functions. 

The system is subject to body forces b(x), heat source r(x) and mass source s(x) which are 
assumed to be C —periodic with period C = [0, L\ x [0,<5L] and to have vanishing mean values on 
C. Since L is a large multiple of e, then C can be assumed to be a representative portion of the 
overall body. This means that the body forces, heat sources and mass sources are characterized by 
a period much greater than the microstructura l size e. 

Following the procedure reported in lBacigalupol ( 2014h . a non-dimensional unit cell Q = [0,1] x 
[0, <5] that reproduces the periodic microstructure by rescaling with the small parameter e is in¬ 
troduced. Two distinct scales are represented by the macrosco pic ( slow) va riables x £ A and 
the microscopic (fast) variable £ = x/e £ Q (see for ex a mple iBakhvalov and Panasenkol ( 1984lf : 
Smvshlvaev and Cherednichenko (2000) and Bacigalupol ( 2014 )1. The constitutive tensors ©, 


m and © are functions of the microscopic variable, whereas the body forces, heat sources and 
mass sources depend by the slow macroscopic variable. Consequently, the mapping of both the 
elasticity and thermodiffusive tensors may be defined on Q as follows: (C^ m ’ £ \x) = C m (£ = 
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as/e), K^ e \x) = K m {£ = x/e), a.^ e \x) = a m {$, = x/e), D^ m ’ e \x) = D m (£ = x/e), /3 {m ’ e \x) 
/3 m (£ = x/e), respectively. 

The relevant micro-fields are the micro-displacement u(x), the microscopic temperature 9(x) = 
T(x)—Tq (To stands for the temperature of the natural state) and the microscopic chemical potential 
rj{x). The micro-stress cr(x), the microscopic heat and mass fluxes q(x) and j(x) are defined by 
the following constitutive relations: 


<r(x) = C m (*) e(x) - (*) 9{x) - f3 m (*) q{x), 

q(x) = -K m (*) V9(x), j(x) = - D m (*) Vq(x), 


(4) 


(5) 


where s(x) = symVit(a;) is the micro-strain tensor which is assumed to be zero at the fundamental 
state of the system. 

Note that, in eqs. ([5]) describing the heat and mass fluxes, we confine ourselves to the essential 
effects and neglect coupling terms, which is an ass umption generally accepted in the quasi-static 
theory of thermodiffusion, see for instance Nowackil 11 974 1. 

The micro-stresses fl]) and the microscopic fluxes ([5]) satisfy the local balance equations on the 
domain A 

V • cr(x) + b(x) = 0 , V • q(x) — r(x) = 0, V • j(x) — s(x) = 0. (6) 

Substituting expressions (Hj)-([5j in equations ([6| and remembering the symmetry of the elasticity 
tensor, the resulting set of partial differential equations is written in the form 

V • ((C m (*) Vu(x)) - V • (a m (*) 0(*)) - V • (/3 m (*) rj(x)) + 6(a) = 0 (7) 

V- (iif m (*) V0(a)) + r(*) = 0, V • (d™ Vtj(®)) + s(a) = 0. (8) 

Moreover, at the interface E between two different phase of the material, the microscopic fields 
satisfy the following interface conditions: 


Ju(a)]] \x££ — 0; 


(C”(j) Vu(a 

= 0 , 


[[0(a)]]| a , gE 

lfo(*)I]|®e£ 


= 0, 


K 

O' 


09W-/3” (?)„(*)) 

= 0 , 


xeY, 


= 0 , 


(?) vfl < x) ■ 
(?) Vv(») • 






= 0, 


(9) 

( 10 ) 

( 11 ) 


where the notation [[/]] = /*(£) — /-' (E) denotes the difference between the values of a function / 
at the interface E separating the phase i from the phase j. 

The micro-displacement, microscopic temperature and chemical potential may be seen in the 
form u(x, £ = x/e), 9(x, £ = x/e), rj(x, £ = x/e) as functions of both the slow and the fast variable. 

It is important to note that since u(x,£), 6(x,£) and 9(x,£) are assumed to be Q— periodic 
smoothing functions with respect to the va riable x . the interfa ce con ditions m-m can be expressed 
directly in function of the fast variable £ ( Bakhvalov and Pa.na.scnkfl 119841 )7 

The solution of microscopic field equations (|7 ]l . ([8 ]) is computationally very expensive and pro¬ 
vides too detailed results to be of practical use, so that it is convenient to replace the heterogeneous 
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model with an equivalent homogeneous one to obtain equations whose coefficients are not rapidly 
oscillating while their solutions are close to those of the original equations. 

Further in the paper, assuming that the size of the microstructure e is sufficiently small with 
respect to the structural size L , an equivalent classical first order thermodiffusive continuum is 
considered. The overall elastic moduli, thermal and diffusion expansion tensors, thermal and diffu¬ 
sive conduction tensors of a homogeneous continuum equivalent to periodic heterogeneous material 
reported in Fig. [T] are derived by means of asymptotic homogenization techniques based on the 
generalization of down-scaling relations. The overall elastic and thermodiffusive properties of the 
homogeneous continuum are expressed in terms of geometrical, mechanical, thermal and diffusive 
properties of the microstructure by means of an asymptotic expansion for the microscopic fields. 
The asymptotic expansion is performed in terms of the parameter e that keeps the dependence on 
the slow variable x separate from the fast one £ = x/e such that two distinct scales are represented. 

In the equivalent homogenized continuum, the macro-displacement U{x ) of component Ui, the 
macroscopic temperature 0(x) and chemical potential T{x) are defined at a point x in the reference 
(e*, i = 1, 2). The displacement gradient is given by VU(x) = = H(x ),and 

then the macroscopic strain is E{x) = symVt/ (x). The macro-stress E(x) associate to E(x) are 
defined as E{x) = Uijei <8> ej with Uij = Sj ,, and the macroscopic heat and mass fluxes are 
respectively: Q(x) = Qiei and J{x) = J^. 


3 Multiscale analysis and asymptotic solution of the hetero¬ 
geneous problem 

3.1 Down-scaling and up-scaling relations 

Follow i ng the approaches developed in [Ba khvalov a nd Panasenkoj ( l984h : ISmvshlvaev and Cherednichenko 
I 2f)0(t) : Bacigaliipo and Gambarottal (12014 1 and lBacigalupo ~ <j2014l l for purely elastic problems in 
periodic heterogeneous media, the microscopic displacement, temperature and chemical potential 
fields are represented through an asymptotic expansion with respect to the parameter e, whose 
terms depend on macroscopic fields and perturbation functions: 


Uk 


(•■«-!) 


+oo 


U k (x) + Y,e l Y, N S q & 


d l U p (a 


1=1 \q\=l 


dx 0 


+oo 

+ E e ' E 

Z=1 \q\=l-l 




dx 0 


= U k (x)+e(Nl 1 J(Z) 


dU p (x) 

dXn-\ 


+ N£\t)Q{x) + N£\z) T{x) 


£=X/e 

£=X/s 


+ e 2 IN' 


+ k: (« 


9 2 u P (x) - r( 2 )f ^de(x) <T( 2) f ^arc{x) 


dXn 


/ £=£C/e 


( 12 ) 
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(-«-!) 


+oo 


t^+E^E K, (0 (o 


2=1 l«l=2 


<9 ; T(a 


fe 0 


£=*/e 






9i / £=CC/e 


w™ (£)EIM 

vv qiq2\^Jf) r 

UX qi UXq 2 / £_ x / £ 


(14) 


In equations (1121) . (1131) and (1141) (commonly known as down-scaling relations), q = qi, 


,qi 


is a multi-index and d l {-)/dx q = d l (-)/dx qi ■ ■ ■ dx qi . Due to their dependence on the slow space 


and W, 


(0 


variable x , the macroscopic fields £4,0 and T are £—periodic functions. Afg , A/g 

are the mechanical, thermal and diffusive fluctuation functions, respectively, whereas N^j and ivg 
denote the additional fluctuation functions corresponding to the contribution of the thermodiffusion 
to local displacement. All these perturbation functions depend on the fast space variable £ = 
x/e, and moreover, a s it will be shown in Section 13.21 they are Q— periodic. Simi larly to the 
procedure reported in Smvshlvaev and Cherednichenkol ( 200(1) and Bacigalupol ( 2014h . the mean 
value of the fluctuation functions is assumed to vanish on the unit cell Q, this means that the 
following normalization conditions are satisfied: 



<im = o, (n$) = \J 

( M $ ) ) = \ j Q M«\z)dZ = 0, 


(t)dZ = 0, (JVg ) = \f (m = 0, 

J w ®{£) d t = 0. (15) 


Introducing a new variable ( £ Q and a vector £ A, which represents the translations of 
the medium with re spect to the £—periodic body forces b(x), heat sources r(x) and mass sources 
s(x) (Bacigal uno and Gambarottal 2014P ). it can be shown that any Q —periodic function g(£ + £) 
satisfies the following invariance property: 


(ff(£ + ^ = \ j Q 9 ^ + ^ di ’ = ^ 9{€ + C)d£. (16) 

According to the invariance property (ll(il) and to the normalization conditions (1151) . the macroscopic 
fields can be defined as the mean values of the microscopic quantities (fl2l) , (fl3l) and (fT4l) evaluated 
on the unit cell Q: 

U k (x) = (u k (x, j +C)) , 0(®) = (^(*. 7 +C)), T(a;) = (rj (x, j +<)) , (17) 

Expressions ED are commonly known as up-scaling relations. More details regarding the structure 
of the down-scaling relations (1121) , (THU) and (THl) are provided in Appendix [D] 
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3.2 First-order asymptotic analysis and derivation of the corresponding 
first-order cell problems 

In order to derive exact expressions for the fluctuation functions affecting the behavior of the micro¬ 
scopic fields Uk,0,il, the down-scaling relations d, Cl and m are substituted into the micro¬ 
scopic field equations 0,0. Remembering the property = j) = (j$f- + = 

( JlT + —) , equation 0 become to the first order approximation 

\ u i / P—cr, / f. 


.-1 


+ Cf, 


(chki* ( ki ), - pi 


Hpq 1 ( x ) + 


T(x) + 




+ bi(x) = 0, 


©O) 

* = 1 , 2 , 


(18) 


where H pqi = dU p /dx qi are the components of the macroscopic displacement gradient tensor pre¬ 
viously defined. Equations 0 assume the following form 


M (1) ■) 
91 j y 

+ 

,* 

K e 

iq i 

de 

dx qi 

H- 

-1- r(x) = 

9i, 3 ^ 

+ 

,i 

D £ 

*9i ,* 

dr 

dx qi 

+ ••• 

-1- s(x) = 


(19) 


( 20 ) 


In order to transform the field equation ©, m and (El in a PDEs system with constant 
coefficients, in which the unknowns are the macroscopic quantities t4(x), 0(x) and T(x), the 
fluctuation functions have to satisfy non-homogeneous equations ( first-order cell problems) reported 
below. 

At the order e _1 from the equation (fT8l) we derive: 


(c-juKIj) 


I fie — 

~ ijpqi,j *P9i’ 


L , (cf jkl N$) . - = nf \ 

( 2 !) 

whereas from thermodiffusion equations Cl and (1201) we obtain: 


\ qi,3) 


+ K e ■ = to (1) (D e -W < ' 1) \ + D E ■ = w (1) 

^ ni qi ’ y^ij VV qi,j) ^ ^iqi,i W qi ’ 


( 22 ) 


where: 


r (1) 

*P9l 


= {Cl 


= 0 , 


= ~( a ij,j) = 0, 


^z (1) = 


m 


£ = (K lA ) = o, w^ = (Dt ai .) = 0. 




9i 


*9i,*/ 


(23) 


The properties (l23ll are consequence of the Q —periodicity of the components Cf,j pqi . a -, , Kf qi 

and Df qi . Note that in equations (fTSl) (E51) the derivatives should be understood in the generalized 
sense. 

The perturbation functions characterizing the down-scaling relations d, d, and d are 
obtained by the solution of the previously defined cells problems, derived by imposing the normal¬ 
ization conditions (11511 . 
















4 Homogenized thermodiffusive Cauchy continuum: field equa¬ 
tions and overall properties 


The field equations of the first order homogeneous continuum can be obtained by the zero order 
terms (equations (IG3D and (KnUl f of the sequence of PDEs derived applying the asymptotic analysis 
to the averaged field equation, see Appendix [A] This implies that the macroscopic displacement, 
temperature and chemical potential are approximated as follows: 

f/ p (x)«t/W(x), e(x)«e<°>(x), T(x)«T‘°>(x). (24) 


Alternatively, the field equations of the equivalent Cauchy continuum can be derived considering 
only the terms of order e° in the equations (15U1) . (RJUl) and (fCTTT) . 

The field equations of an homogeneous first order continuum in presence of thermodiffusion are 
given by 


a 


9 2 Ur, 


™ dx qi dx q2 


x iqi 


<90 

dx, 




it 3l 


9l 


dr 

dx 


+ bi — 0, 


9l 


(25) 


K, 


<9 2 0 


9192 dx qi dx q2 


+ r = 0, D, 


d 2 T 


9192 dx qi dx q2 


+ 5 — 0, 


(26) 


where Ci qipq2 are the components of the overall elastic tensor, ai qi and /3i qi are respectively the 
components of the overall thermal dilatation and diffusive expansion tensors, K qiq2 denotes the 
components of the overall heat conduction tensor and D qiq2 represents the components of the overall 
mass diffusion tensor. Remembering the approximation (12411 . the macroscopic field equations (1251) - 
(1251) can be compared to the zero order terms of the averaged field equation (16311 and (1M1) for 
determining the overall properties of the thermodiffusive Cauchy continuum. In order to relate the 
coefficients , h,-^, nj q j , m^ q2 , w q l\ 2 contained in the equations (l63l) and (l64l) to the overall 

elastic and thermodiffusive constants of the media Ci qipq2 , ai qi , pi qi , K qiq , 2 , D qiq2 , the symmetries 
of the tensors of components n\ ' q2 , n iqi , fi iqi , m qi q2 , w qi q2 , and the ellipticity of the field equations 
(l63l) and (l64l) are required. A demonstration of these properties is reported in Appendix ICl As a 
consequence of these properties, it can be observed that: n^ qiq3 = \(Ci qi pq 2 + Ciq 2 pq \), h[ q ^ = cti qi , 
™iql = Piqn m 9 i 92 = Rq l92 and w’gik = ^ 9192 - I 11 particular, comparing the field equation (E5l) to 


and remembering the relationship between n\ vq 

d 2 U 

repetition of the indexes <71 and 52 : Cit 


z pqi 92 an< ^ ^*9iP92 < ^ is easy to note that due to the 


y &9iP92 dx qi dxq 2 n ip9i92 dx qi dx q2 ~ 2 V^*9iP92 ~v'-^92P9i 1 dx n dx q2 ' 

The overall elastic and thermodiffusive tensors, obtained in terms of fluctuation functions, and 
the components of microscopic elastic and thermodiffusive tensors, take the form (see Appendix [A] 
for details): 


( 2 ) 


d'Ur, 


= h(c i( 


> rCi, 


<); 


d z U v 


Ci qiP q 2 — {C r j kl 


{ C rjkl ( 


N, 


(i) 

riqi ,j 


idjqi 


■N 


(i) 


rqn,J 


S rqi dij 


<L 




+ $kq2 &lp 


)) 



The components Ci qipq2 , K qiq2 and D qiq2 of the overall constitutive tensors of the material co¬ 
incide with those derived by asymptotic homogenization techniques applied to uncoupled static 
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elastic ( Bakhvalov and Panasenkol 1984: Smvshl vaev and Cherednichenkol 2000; Bac igalupol 2014 1 
and heat conduction problems ( Zhang et all l2007t) in media with periodic microstructures. The 
components on qi and f3i qi of the coupling thermodiffusive tensors have been obtained by means 
of a consistent generalization of the down-scaling relations (fT2l) (fl3l) and (ED. These expressions 
relate the microscopic displacement field to the macroscopic displacements, temperature, chemical 
potential and their higher order gradients. 


5 Illustrative example: homogenization of bi-phase orthotropic 
layered materials in presence of thermodiffusion 

The general results obtained are now applied to the case of a bi-phase layered material in presence of 
thermodiffusion. Exact analytical expressions for the overall elastic and thermodiffusive constants 
are derived. Considering a two-dimensional infinite thermodiffusive medium subject to periodic 
body forces, heat and/or mass sources, the solution obtained applying the proposed homogenized 
model is compared with the results provided by the analysis of the corresponding heterogeneous 
problem. 




Figure 2: (a) Heterogeneous and homogenized models with C -periodic body force bj, heat sources 
'r(xj) and mass sources s(xj); (b) Periodic cell and constituents: bi-phase layered material. 


5.1 Perturbation functions and overall constitutive constants: exact an¬ 
alytical expressions 

Let us consider a layered body obtained as an unbounded cfe— periodic arrangement of two different 
layers having thickness a and b , where d ,2 = s = a + b and </ = a/b are defined. The phases 
are assumed homogeneous and orthotropic, with an orthotropic axis coincident with the layering 
direction ei, the geometry of the system is shown in Fig. [2] The orthotropic symmetry is assumed 
for both the elastic and thermodiffusive tensors. The micro-fluctuation functions N^ qi , N^\ 

M q l^ and W q ^ are analytically obtained through the solution of the cell problems formulated in 
Section 13.21 (see equations (f2TTl and (l22l) and conditions (l23l) ). Due to the particular properties of 
symmetry of the microstructure, these functions depend only on the fast variable £2 ■ This variable 
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is perpendicular to the layering direction ei (see Fig. 0). The non-vanishing micro-fluctuation 
functions and N^\ obtained by solving the cell problem of order e - 1 m are: 


AT {l)_ a _ ^1122 °1122 ea . 

iV 2ii ~ ra _ r rt> ^ 2 ’ 

°2222 W2222 

. . r^b s~hl 
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N^\- b = C 1122 “ 22 & 
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(28) 


N^- a = 


N^~ a = 


x 22 


x 22 


0 2222 + CC*22 


-C 2 a ; ^V 2 (1) - b = C 


P 22 ~, 


C 2222 + CC 2: 


2222 
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x 22 
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(29) 

(30) 


where £ 2 € 


< C 
2(C+1)’ 2(C+1) 


and £2 € 


1 1 

2(C+1) ’ 2(C+1) 


are non-dimensional vertical coordinates 


centered in each layer. The non-vanishing fluctuation functions associate with the thermodiffusion 
equations, derived by the solution of the cell problems of order £ _1 (1^1) are: 


M 2 (1) - q = - 


K22 - I<22 ,q 

iq 2 -a<!k 


M (D_ b = iq 2 K 22 g 
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^22 C-O 


(32) 


Note that the superscripts a,b denote that the elastic and thermodiffusive constants are referred 
respectively to the phases a and b. 

In order to derive the overall elastic and thermodiffusive constants corresponding to a first 
order equivalent continuum, the fluctuation functions dsn), m , (ra, m and (f32|) are used into 


expressions Ezl). The components of the overall elastic tensor Ci qivq2 take the form: 
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( 33 ) 

The non-vanishing components of the thermal dilatation tensor ai qi and diffusive expansion tensor 
(3i qi are respectively given by 


an = 


C(Oll22 a 22 ~ Oll22 a 22 ~ Q2222 a ll Oll22 a 22 + Oll22 a 22 ^2222 a ll) ~ C~C-2222 a ll O2222 a ll 
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(34) 
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The non-vanishing components of the heat conduction tensor K lqi and mass diffusion tensor Di qi 


take the form 


K-n = 
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CK 11 


C + l 


D„ = 
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11 


(D 


11 


C + l 


A 22 = 


D 22 = 


(C+l )K&K, 


K%2 + c+1 2 : 
(C + 1)^12^22 


n« 


C^2 


(36) 

(37) 


Considering the case of isotropic phases, the components of the elasticity tensor become Cl lxl = 
rv ' — ZsMs. C{ 212 = 2 (Boa ’ (^H 1 ? = a, 6), where for plane-strain: E q = 


C2222 


C1122 


1-0 f’ w 1122 l-z/^> '“1212 2(1+S C )’ “’'Vi I""**"- 1 - 

v q = , whereas for plane-stress: E q = v q = i/ ? , being E q the Young’s modulus and 

the Poisson’s ratio, respectively. The components of the thermal dilatation and diffusive expansion 
tensors take respectively the forms: a\ x = a 22 = cP, 0^ = 022 = 0““ (note that the coefficients cP 
and 0 ? can be expressed in te rms of the li near isotropic thermal and diffusive expansion coefficients 
and the elastic moduli ( Nowacki . 1974 . 1986l l. The components of the heat conduction and mass 


diffusion tensors finally become A')-, = K 22 = AT 5 and D\ ± = D 22 = £>+ The overall elastic and 
thermodiffusive constants for the case of isotropic phases are reported in Appendix [C] 

By an asymptotic expansion of the constants (1551) . (1551) . (1551) . (155)l and (1571) in terms of the 
concentration of the two constituents phases (not reported here for conciseness), it can be easily 
shown that, if the concentration of the phase a vanishes, the overall elastic and thermodiffusive 
constants of the bi-phase layered material tend to the values corresponding to phase b. Conversely, 
if the concentration of the phase a tends to one, the same expressions tend to the elastic and 
thermodiffusive constants of the phase a. 

In order to simplify the required computations, for the illustrative examples both the phases 
are assumed to be isotropic, and then the overall elastic and thermodiffusive constants reported in 
Appendix [C] are used. These constants can be represented in the non-dimensional form: 


^7»9lP92 (pC ) C) T'a j Vb ) — 
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(38) 


'■iqi ^iqi 

lb _ pa /pb _ _ iPa ! T>b „„ _ n« / f>& 


where p c = E a /E bl p a = a a /a b : pp = 0 a /0 = K a /K b ,p D = D a /D b , and C iqipq2 = (Cf qipq3 + 

C b qipq2 )/2,a iqi = « 91 +a b qi )/2j iqi = (0“ gi + f3 b qi )/2, K iqi = {Kf qi + AT^J/2, D iqi = (D“ qi + 
D b )/2. It is important to note that if the Poisson’s coefficients of the two phases are identical 


(i.e. v a = i/;,), the non-dimensional overall elastic and thermodiffusive constants 
following property: 


possess the 


C‘iq 1 pq 2 (PCi Cj V) — ^iqipq 2 (Pc )C i ^) > &iqi(PC ) Pen Ci &) — ai qi (pp , P a , C >^)> 


Piqiipc 1 Pen C>v) — Pi qi (pc ) P a jC > ^)i &iqi (PK: C) — ATi 9l (p K , C )j 
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( 39 ) 


Di qi {pa, £) — D iqi (p K >C )■ 

The variation of the normalized components of the overall elasticity tensor Cnn and C2222 with 
the ratio pc is reported in Figs. [3j/(a) and|4]/(a), respectively. The same value of the Poisson’s 
coefficient v = 0.3 has been assumed for both the phases, and several values of the non-dimensional 
geometrical parameter £ has been considered for the computations. It can be observed that for 
pc = 1, corresponding to the case of two isotropic phases having identical elastic properties, the non- 
dimensional components of the overall elastic tensor assume the value Ci qipq2 = 1 (i. e. Ci qipq2 = 
Ci qiP q 2 = (-'i qiP q 2 )- In Figs. [3j/ (6) and0F(6) the components Cnn and C2222 are plotted as 
functions of £ for different values of pc considering the fixed Poisson’s coefficient v = 0.3 identical 
for both the phases. For £ —> 0, the thickness of the phase a vanishes. Consequently, the values 
of the overall elastic constants tends to those of the phase b: Ci qipq2 = Cf qipq2 , and the limit 
values assumed by the normalized components of the elastic tensor reported in the figures are 
Ciq ipq2 = 2/(1 + pc)- Conversely, for £ —> +00 the thickness of the phase b tends to zero, and then 
Ci qi pq 2 = Ci qipq2 and the non-dimensional constants plotted in Figs. EF(fr)-SF W assume the limit 
values Ci qipq2 = 2pc/iX + pc)- Results (not reported here) show that the normalized components 
C 1212 and C 1122 have a behaviour very similar to that of the component 02222 - 



Figure 3: (a) Dimensionless constant Cun vs. the ratio pc for v a = f'h = 0.3 and for different 
values of the geometric ratio £: £ = 1/10 green line, £ = 1/5 blue line, £ = 1/2 red line, £ = f/f 
black line, £ = 2 red points, £ = 5 blue points, £ = 10 green points, (b) Dimensionless constant 
Cnn vs. the geometric ratio £ for v a = Ub = 0.3 and for different values of the ratio pc- pc = 2 
red line, pc = 5 blue line, pc = 10 green line, pc = 30 black line, pc = 50 violet line. 
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Figure 4: (a) Dimensionless constant C2222 vs. the ratio pc for i' a = = 0.3 and for different 

values of the geometric ratio £: £ = 1/10 green line, £ = 1/5 blue line, £ = 1/2 red line, £ = 1/1 
black line, £ = 2 red points, £ = 5 blue points, £ = 10 green points, (b) Dimensionless constant 
C2222 vs. the geometric ratio £ for v a = i/f, = 0.3 and for different values of the ratio pc' pc = 2 
red line, pc = 5 blue line, pc = 10 green line, pc = 30 black line, pc = 50 violet line. 


The three-dimensional plots reported in Fig.[5]show the variation of the normalized components 
of the overall thermal dilatation tensor an and 0:22 as functions of p c and p a , assuming v = 0.3 
for both the phases and £ = 1. In Figs. EJ/(a) and[7j/(a) the variation of dn and 0:22 with the 
non-dimensional ratio pc is reported for several values of £ assuming v a = ui, = v = 0.3 and p a = 2. 
For pc = 1, corresponding to the case of two isotropic phases with identical elastic constants but 
different thermal dilatation properties, the normalized components of the overall thermal dilatation 
tensor tend to the values a iqi = 2((p a + 1 )/[{p a + 1)(C + 1)] (be. a igi = « ?1 £ + a- ?1 )/(C + !))• 
In the case where pc = 1 and also p a = 1, both the elastic and thermal dilatation tensors of the 
two phases are identical, and then &i qi = 1. In Figs. and [7 \{b) the same constants dn 

and &22 are plotted as functions of £ for u = 0.3, p a = 2 and several different values of pc- For 
£ —» 0, the thickness for the phase a vanishes, and the elements of the overall thermal dilatation 
tensor tends to those of the phase b (i.e. ai qi = a\ qi )■ As it can be observed in the figures, in 
this case the normalized constants tend to a limit value which is the same for any value of pc (i.e. 
d i qi = 2/(1 + p a )). This value can be easily derived by using expressions for an and <222 reported 
in Appendix [C] Conversely, for £ —> + 00 , the thickness of the layer b tends to zero, the effective 
thermal dilatation constants tend to those of the phase a (i.e. ai qi = o;“ ) and the normalized 
components dn and 0122 reported in Figs. [6j/(6) and[7]/(6) assume the values d; 9l = 2p a /{^ + p a )- 
The properties of the normalized elements of the overall diffusive expansion tensor /3n and P 22 are 
similar to those of dn and 0.221 an d can be easily studied substituting the non-dimensional ratio 
Pa with pp . 

The variation of the normalized components of the overall heat conduction tensor and 
K 22 with the non-dimensional ratio px are shown in Figs. |S]/(a) and[!2J/(a). Several values of £ 
have been assumed for the computations. It can be observed that for px = 1, we have Ki qi = 1. 
This is due to the fact that the value px = 1 corresponds to the case where the heat conduction 
of the two phases are identical, and then Kf = K\ . In Figs. [8j/ (a) and [9]/ (fo) the same non- 
dimensional components A'n and K 22 are reported as functions of £ for different values of px ■ As 
it is shown by these figures, for £ —»• 0, I\n and K 22 tends to a finite value depending on px- This 
limit correspond to the case of vanishing thickness of the layer a, where Ki qi = K ^ qi , and then 
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Figure 5: Dimensionless component an (a), 0:22 (b) vs. the ratios pc and p a for u a = Vb = 0.3 and 

C = 1/2. 



Figure 6: Dimensionless component an vs. the ratio pc for v a = Vb = 0.3, p a = 2 and for different 
values of the geometric ratio £: £ = 1/10 green line, £ = 1/5 blue line, ( = 1/2 red line, ( = 1/1 
black line, ( = 2 red points, ( = 5 blue points, ( = 10 green points, (b) Dimensionless component 
an vs. the ratio pc for v a = f'h = 0.3, p a = 2 and for different values of the geometric ratio pc'- 
pc = 2 red line, pc = 5 blue line, pc = 10 green line, pc = 30 black line, pc = 50 violet line. 


Kiq 1 = 2/(1 + pk )■ In the limit / —»• + 00 , for which the thickness of the layer b vanishes and 
Ki qi = Kf , the normalized components of the overall heat conductivity tensor tend to a finite 
value given by Ki qi = 2px/{]-+ Pk)- The non-dimensional components of the overall mass diffusion 
tensor Dn and D 22 are characterized by properties similar to those of Kn and K 22 , and can be 
studied substituting the non-dimensional ratio pk with pc- 
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Figure 7 : Dimensionless component <222 vs. the ratio pc for v a = Vb = 0.3, p a = 2 and for different 
values of the geometric ratio £: £ = 1/10 green line, £ = 1/5 blue line, ( = 1/2 red line, ( = 1/1 
black line, £ = 2 red points, £ = 5 blue points, ( = 10 green points, (b) Dimensionless component 
(222 vs. the ratio pc for v a = Vb = 0.3, p a = 2 and for different values of the geometric ratio pc' 
pc — 2 red line, pc = 5 blue line, pc = 10 green line, pc = 30 black line, pc = 50 violet line. 



Figure 8: Dimensionless component A '11 vs. the ratio px for different values of the geometric ratio 
(: C = 1/10 green line, ( = 1/5 blue line, ( = 1/2 red line, £ = 1/1 black line, ( = 2 red points, 
( = 5 blue points, ( = 10 green points, (b) Dimensionless component K\\ vs. the ratio px for 
different values of the ratio px'- Pk = 2 red line, px = 5 blue line, px = 10 green line, px = 30 
black line, px = 50 violet line. 
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Figure 9: Dimensionless component K 22 vs. the ratio px for different values of the geometric ratio 
(: £ = 1/10 green line, £ = 1/5 blue line, £ = 1/2 red line, £ = 1/1 black line, £ = 2 red points, 
£ = 5 blue points, £ = 10 green points, (b) Dimensionless component K 22 vs. the ratio px for 
different values of the ratio px- px = 2 red line, px = 5 blue line, px = 10 green line, px = 30 
black line, px = 50 violet line. 
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5.2 Comparative analysis: homogenized model vs heterogeneous mate¬ 
rial 


In order to study the capabilities of the proposed homogenization procedure, the two-dimensional 
bi-phase orthotropic layered material shown in Fig. [2] is assumed to be subjected to ^-periodic 
harmonic body forces hi, directed along the orthotropy direction Xj (see Fig. ED and Zl-periodic heat 
and mass sources r(xj) and s(xj ): 


2nmxj 2iniXj 2 npxj 

i - ?- % - 

bj(xj) = Bje Lj ( r(xj) = Re ? s(xj) = Se ; (40) 

where: j = 1, 2; L\ = L 2 = L\ Bi, R,S £ C; m,n,p € Z; and i 2 = — 1. 

This problem is analyzed by applying the homogenized first-order model with overall elastic and 
thermodiffusive constants derived from the homogenization of the periodic cell through the approach 
developed in previous Sections. The obtained results are then compared with those derived by means 
of a fully heterogeneous modelling procedure. Due to the periodicity of the heterogeneous material, 
body forces, heat and mass sources considered, only an horizontal (or vertical) characteristic portion 
of length L of the heterogeneous model is analyzed (Fig. [ 2 ) 3 ). In order to assess the reliability of 
the homogenized model, the macroscopic displacement, temperature and chemical potential fields 
are compared to the corresponding fields in the heterogeneous model by means of the up-scaling 

relations ed- 

The overall elastic and thermodiffusive constants involving the fluctuation functions are obtained 
in exact analytical forms via expressions ED, ED, (ED, ED, and (1361) . Conversely, the solution 
of the heterogeneous problem with £—periodic harmonic body forces is computed via FE analysis 
with periodic boundary conditions on the displacement temperature and chemical potential fields. 
For the considered two-dimensional body subject to body forces along the orthotropy axes, heat 
and mass sources, the homogenized field equations (1251) - (1261) take the form: 


C; 


d 2 U, 

' 9x 2 
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— Qtj 


-h 
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n dxi n dxi 


+ bn = 0, K 


5 2 0 


n dx 2 


+ r = 0, D, 


d 2 r 


s = 0, 


( 41 ) 


where j = 1,2 are not summed indexes. Equations (1411) describe an extensional problem in pres¬ 
ence of thermodiffusion. Considering body forces, heat and mass sources of the form (l40ll . the 
macroscopic displacements, temperature and chemical potential fields are given by 


Uj(xj) = 


BjL 2 


Cjjjj {2irm) 2 


0 —l 


RotjjL 3 
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( 43 ) 


where j = 1,2 are still not summed indexes. In order to compare the behavior of the derived 
analytical solution with the numerical results provided by the heterogeneous model, only the real 
part of macroscopic fields (ED, (ED accounted. Moreover, the imaginary part of the amplitudes 
Bj , R and S is assumed to be zero. The real part of expressions (1421) . (H51) can be written in the 
non-dimensional form: 
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(44) 
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Q*(xj) = cos 


(2nnXn 


Lj 


T *(xj) = cos 


( 2 npXn 

{— 


where U* = 


Uj Cjjjj ( 27 rm) 
B~L 3 
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0A' JJ (27rw) 2 
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rD J3 (2Kp) 2 
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and ££ = 


a jj RLj 

Knn Bn 


and s£. = 


(45) 


Pjj SLj 

DnnBn ’ 


j = 1,2 are still not summed indexes. 

The amplitude functions E and E^ are associated respectively with the thermal expansion 
and mass diffusion contribution to the macroscopic displacement (14411 along the direction e 7 . In 
order to study the influence of the geometrical, elastic and thermodiffusive properties of the phases 
on E “■ and E^-, the following non-dimensional form for these functions is introduced: 


£7 cm 

S fj(PC,Pa,PK,C,Pa,Pb) = -=f, -jj{PC, P/3, PK,C.,Pa,Pb) = with j = 1, 2, (46) 

“li ^jj 

where E “ = RLj (a a + a b )/ Bj (. K a + JtT b ) and E^ = RLj (/?“ + /3 fe )/(Z? a + D b ). In the case where 
the two phases possess the same value of the Poisson’s coefficient ( v a = Pb ), the following property 
is verified for the normalized amplitude functions: 

jj^PC ’ Pat PK t Cl ^) 1 — jj(Pc ’Pa ’Pk l —‘ jj (PC : P3 i PK ; C: ^0 1 — jj^Pc ’Pa ’Pk (1^) 

The three-dimensional plot reported in Fig. HO^fq) shows the variation of the normalized ampli¬ 
tude component E with p a and px for P a = Pf, = 0.3, pc = 10 and £ = 1. In Fig. I10[/fbl the 
same component , which represents the contribution of the thermal expansion to the macro¬ 
scopic displacement along the direction ei, is plotted as a function of pc assuming P a = Pb = 0.3, 
Pa = Pk = 10 and several values of the dimensionless ratio C ). The variation of the normalized 
component E 22 , corresponding to the contribution of the thermal expansion to the macroscopic 
displacement along e 2 (see Fig. [2]), is reported in Fig. W(°) as a function of p a and px for 
P a = Pb = 0.3, pc = 10 and C, = 1, and in Fig. ITT)/ f61 as a function of pc assuming P a = Pb = 0.3, 
Pa = Pk = 10. Observing the curves reported in the figures, it can be noted that the normalized 
amplitude E 22 , associate to the component of the macroscopic displacement U 2 {x 2 ) parallel to 
the stratification direction, is greater than the amplitude E^ which correspond to the component 
Ui(x i) parallel to the stratification direction. The dimensionless amplitude and E 22 , associ¬ 
ated to the contributions of the mass diffusion respectively to U^(x i) and U 2 {X 2 ), are characterized 
by the same properties of H’“ 1 and E 22 , and their behavior can be easily studied substituting the 
non-dimensional ratios p a and px with pp and pc- 
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Figure 10: (a) Dimensionless amplitude S“ x vs. the ratios p a and px for v a = = 0.3, pc =10 

and £ = 1. (b) Dimensionless amplitude vs. the ratios pc for v a = Vb = 0.3, p a = px = 10 and 
£ = 1 for different values of the geometric ratio /: £ = 1/10 red line, £ = 1/5 blue line, / = 1 black 
line, £ = 5 blue line points, £ = 10 red line points. 



Figure 11: (a) Dimensionless amplitude vs. the ratios p a and px for u a = Db = 0.3, pc= 10 
and £ = 1. (b) Dimensionless amplitude E% 2 vs. the ratios pc for v a = f'b = 0.3, p a = pif = 10 and 
£ = 1 for different values of the geometric ratio £: £ = 1/10 red line, £ = 1/5 blue line, £ = 1 black 
line, £ = 5 blue line points, £ = 10 red line points. 
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The analytical solution (Hil) and (H5l) , derived by the solution of the homogenized field equations 
EID is now compared with the results obtained by the finite element analysis of the heterogeneous 
problem corresponding to the bi-phase layered material reported in Fig. [2] subject to harmonic body 
forces, heat and mass sources. More precisely, finite element analysis of the heterogeneous problem, 
performed by means of the program COMSOL Multiphysics, provides the local fields Uj, 6 , rj which 
are used together with the up-scaling relations (ED for obtaining the macro-scopic fields Uj, 0 
and T. These macro-scopic quantities are compared with the analytical expressions (HU) and fi5l) . 
Plane stress condition has been assumed for both the solution of the homogenized equations and 
the heterogeneous problem, and two isotropic phases with the same value of the Poisson’s coefficient 
u a = i*b — 0.3 have been considered. 

In Fig. [HI the macroscopic displacement component Uf and temperature 0 * evaluated using 
analytical expressions (HU) and fl5l) \ are reported as functions of the normalized spatial coordinate 
x\/L (continuous lines in the figure) and compared with the numerical results obtained by the 
heterogeneous model assuming periodic body forces b±(xi) and heat sources r{x i) and considering 
the value of the amplitude S{\ = 1 (diamonds in the figure). The following values for the geometrical 
parameters, the ratios between the elastic and of thermodiffusive constants have been assumed: 
L/e = 10, pc = 10, p a = 10, Pi_ 3=0 pk = 10, pd = 0, the effects of the mass diffusion have been 
neglected in this example. The macroscopic displacement and temperature fields are plotted for the 
characteristic portion of length Li = L , corresponding to X\/L = 1 (i. e. for 0 < X\/L < 1), and 
several values for the wave numbers m,n £ Z have been considered. Observing the curves, for both 
the quantities t/f (aq/L) and 0*(aq /L) a good agreement is detected between the results derived by 
means of the first order homogenization approach and those obtained by the heterogeneous model. 

Results for the macroscopic displacement component U$ and temperature in 0* along the char¬ 
acteristic portion of length L 2 = L in direction parallel to e 2 (not reported here for conciseness) 
show a good agreement between the solution obtained by means of the first-order asymptotic homog¬ 
enization method and the values obtained by means of finite element analysis of the heterogeneous 
problem. 




Figure 12: Macroscopic displacement component Ul and temperature fields 0* due to harmonic 
body force b\ and temperature source r. The heterogeneous model (diamonds) is compared with 
the homogenized first order model, (a) Dimensionless macro displacement U* vs. the ratio xi/L 
for = 1 for different values of wave number n, m (n = 1, m = 1 red line; n = 2, m = 1 blue 
line; n = 1, m = 2 green line), (b) Macro temperature 0* vs. the ratio x\/L for different values of 
wave number n (n = 1 red line; n = 2, blue line). 
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Figure 13: Macroscopic displacement component JTj*, temperature 0*, and chemical potential T* 
due to harmonic body force b\ and temperature and mass sources r,s respectively. The heteroge¬ 
neous model (diamonds) is compared with the homogenized first order model. Dimensionless macro 
displacement U* vs. the ratio x\/L for = 1 , 5 ^ = 1 and wave number m = 1 (a), m = 2 (b) 
for different values of wave number n, p (n = 1, p = 1 red line; n = 2, p = 1 blue line; n = 2, 
p = 2 green line), (c) Macro temperature 0* vs. the ratio X\/L for different values of wave number 
n (n = 1 red line; n = 2, blue line), (d) Macro concentration T* vs. the ratio x\/L for different 
values of wave number p (p = 1 red line; p = 2, blue line). 


In Fig. [H the variation of the normalized component of the macroscopic displacement Ui(x\/L), 
temperature Q*(x± /L) and chemical potential T*(xi /L) along the characteristic portion of length 
Li = L is plotted as a function of x\/L. Two isotropic phases having the same Poisson’s coefficient 
v = 0.3 have been assumed, and the same values of the previous example have been assigned to the 
geometrical, elastic and thermal parameters. The amplitude of the mass diffusion contribution to 
the displacement is assumed to be = 1, and px = 10. Similarly to the previous case, for the 
finite element analysis of the heterogeneous elastic thermodiffusive problem harmonic body forces 
6i(xi), heat and mass sources r{x i) and s(a;i) have been introduced. The reported curves show a 
good agreement between the results obtained by the asymptotic homogenization (continuous lines 
the figure) and those provided by the heterogeneous elastic thermodiffusive model (diamonds in 
the figure). The good agreement between the results coming from the two different approaches 
can be observed for U*{x\/L) in Fig. [13]/ (a), (6), for Q*(x\/L) in Fig. [13]/ (6) and for T*(x\/L) in 
Fig.[T5|/fd'). Similar results are obtained for the macroscopic displacement l/j , temperature 0* and 
chemical potential T* along the characteristic portion of length L 2 = L in direction parallel to e 2 
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6 Conclusions 


A general asymptotic homogenization approach for describing the static elastic, thermal and diffu¬ 
sive properties of periodic composite materials in presence of thermodiffusion is proposed. Down- 
scaling relations associating the displacements, temperature and chemical potential at the micro¬ 
scale to the corresponding fields at the macro-scale are introduced. Perturbation functions are 
defined for representing the effects of the microstructures on the microscopic displacement, tem¬ 
perature, chemical potential and on the coupling effects between these fields. These perturbation 
functions are obtained through the solution of non-homogeneous problems on the cell defining 
periodic boundary conditions and normalization conditions ( up-scaling relations). 

Averaged field equations of infinite order are derived for the considered class of periodic thermod- 
iffusive materials, and an original formal solution is performed by means of power series expansion 
of the macroscopic displacements, temperature and chemical potential fields. Field equation for the 
homogenized Cauchy thermodiffusive continuum are derived, and exact expressions for the overall 
elastic and thermodiffusive constants of this equivalent first order medium are obtained. 

An example of application of the developed general method to the illustrative case of a two- 
dimensional bi-phase orthotropic layered material is provided. The effective elastic and thermodif¬ 
fusive constants of this particular composite material are determined using the general expressions 
derived by the asymptotic homogenization procedure. Analytical expressions for the macroscopic 
fields derived by the solution of the homogenized equations corresponding to the first order equiv¬ 
alent continuum. Finite element analysis of the corresponding heterogeneous model is performed 
assuming periodic body forces, heat and mass sources acting on the considered bi-phase layered 
composite. In order to compare the analytical solution of the homogenized equations with the 
numerical results obtained by the heterogeneous model, the microscopic fields computed by finite 
elements techniques are used to estimate the macroscopic displacements, temperature and chemical 
potential fields by means of the up-scaling relations defined in the paper. The good agreement 
detected between the solution derived by the homogenized first order equations and the numer¬ 
ical results obtained by the heterogeneous model through the up-scaling relations represents an 
important validation of the accuracy of the proposed asymptotic homogenization approach. 

Thanks to the great versatility of the asymptotic homogenization techniques and to the pro¬ 
posed general rigorous formulation, the method developed in the paper can be adopted for studying 
effective elastic and thermodiffusive properties of many composite materials, without any other as¬ 
sumption regarding the geometry of the microstructures in addition to the periodicity. In particu¬ 
lar, the proposed asymptotic homogenization procedure can have relevant applications in modelling 
mechanical and thermodiffusive properties of energy devices with layered configurations, such as 
lithium ions batteries and solid oxide fuel cells. In standard operative situations, the components 
of these devices are commonly subject to severe thermomechanical and diffusive stress which can 
cause damages and crack formation compromising their performances. Consequently, evaluating the 
overall elastic and thermodiffusive properties of these battery devices through the asymptotic ho¬ 
mogenization approach illustrated in the paper can represent an important issue in order to predict 
damaging phenomena and to improve the efficient design and manufacturing of these systems. 

Multi-scale homogenization techniques such as that proposed in the paper provide an accu¬ 
rate description of the macroscopical mechanical and thermodiffusive properties of heterogeneous 
materials through the derivation of effective constants of the first order equivalent continuum. Nev¬ 
ertheless, first order homogenization procedures are not enough accurate to model size-effects and 
non-local phenomena connected to the microstructural scale length. As a consequence, the devel- 
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oped first order homogenization approach does not provide a precise description of the behavior of 
thermodiffusive composite materials in presence of high gradients of stresses, deformations, temper¬ 
ature, chemical potential, heat and mass fluxes, as well as of non-local phenomena such as waves 
dispersion. In order to overcome these limits in the accuracy, non-local higher order homogenization 
techniques can be used. These methods provide constitutive relations of equivalent higher order 
continuum media including characteristic scale-lengths associate to the microstructural effects. Us¬ 
ing the rigorous and general approach illustrated in the paper, a better approximation of the elastic 
and thermodiffusive behavior of composite materials in presence of strong gradients can be obtained 
through the solution of higher order cells problems involving the coefficients of the averaged field 
equations of infinite order reported in Appendix [A] As it is shown in the same Appendix, these 
equations can be formally solved by means of a double asymptotic expansion performed in terms 
of the microstructural size. 
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A Higher-order analysis and averaged equations of infinite 
order 

In this Appendix, explicit expressions for the higher order cells problems associated to the down- 
scaling relations (IT21) , (fT51) and (1141) are reported. Moreover, the averaged field equations of infinite 
order are derived, and a formal solution is obtained by means of an asymptotic expansion of the 
macroscopic fields in terms of the microstructural size. 


A.l Higher-order asymptotic analysis and derivation of the correspond¬ 
ing cell problems 

In order to derive exact expressions for the fluctuation functions affecting the behavior of the micro¬ 
scopic fields Uk,0, g, the down-scaling relations (fl^l) . (TT51) and (11411 are substituted into the micro- 
pic field equations 0, ®. Remembering the property g^-/(cc,£ = j) = 
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where H pqi = dU p /dx qi are the components of the macroscopic displacement gradient tensor pre¬ 
viously defined, and n pqiq2 = d 2 U p /dx qi dx q2 are the elements of the macroscopic second gradient 
tensor. Equations (j8]) assume the following form 



+.+ r(x) = 0, 
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In order to transform the field equations (H5|) . C^l) and (15(11) in a PDEs system with constant 
coefficients, in which the unknowns are the macroscopic quantities U k (x), 0(x) and T(x), the 
fluctuation functions have to satisfy non-homogeneous equations ( cell problems ) reported below. 

At the order e^ 1 from the equations m, mi, m we derive the first-order cell problems 
reported in Sec. 13.21 in the text of the paper, equations m and (1221) . 

At the order e , equation (|48l) yields the following second-order cell problems 
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At the same order, from (1491) . (1501) we derive the second-order thermodiffusive cell problems: 
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In general, for the m = 1, 2, • • •, the to— order cell problems associate to equation (1481) assume the 
form: 
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whereas the to— order thermodiffusive cell problems corresponding to equations (1491) and (1501) are: 
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The perturbation functions characterizing the down-scaling relations ©, (USD, and da are 
obtained by the solution of the pre viously defined cells problems, de rived by imposing the normaliza- 
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tion conditions (1T51) . According to Bakhvalov and Panasenko ( 

(2000), the constants (1541) and (1551) are determined by imposing that the non-homogene ous terms in 


equations (1551) . (15I1) . (1571) . (1551) . (1551) and (1571) (associated to the auxiliary body forces ( Bacigalupoi , 
201 4). heat and mass sources) possess vanishing mean values over the unit cell Q. This implies 
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continuity and regularity of the microscopic fields (micro-displacements, micro-temperature and 
micro-concentration) at the interface between adjacent cells are guaranteed. 


A.2 Averaged field equation of infinite order and its formal solution 

Using the cell problems (17T1) . (1771) . (1551) . (1511) . (1571) . (1551) . (1551) and (1571) together with the constants 
definitions (1751) , (1541) and (1551) into the microscopic field equations (fl8l) , (fl9l) and (1201) , the averaged 
equations of infinite order are derived: 
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where q is a multi-index, d n+1 {-)/dx q = d n+1 (-)/dx qi ■ ■ ■dx qn+j with j € N, 
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A formal solution of the averaged field equations of infinite order (l59l) . (IfiOl) and (RTTh is ob¬ 
tained by means of an asymptotic expansion of the macroscopic fields Ui , 0 and T in terms of the 
microstructural size e, i.e. 
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By substituting the series (EH) into (1591) , (1601) and (1611) , a sequence of equations for determining 
the terms of the asymptotic expansion u[ m \ 0( m ) and is obtained. At the order e°, from the 
equation (1551) we derive: 
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At the generic order m from (1591) we obtain 


(63) 


(64) 


3 (2) C/^ m) „ (2 )<90( m ) 

- J_ 77 : ’ - 

ip9192 dx qi dx q2 191 dx, 


A-2) 


,( 2 ) 


9i 


^Y(m) 

dx, 


ra +2 


qi 


EE- 

r=3 |q-|=r 


(r) 9 r Up 


(ra+2—r) 


m+2 


V- V- (r ) 9 r - 1 0( m + 2 - r ) ’^+ 2 „ i r - l a r ' 1 T( m + 2 - r ') 

E E ^— + E E — -= 0 ' 


* 9 

r=3 |g|—r—1 

and (l60l) and (lUTl) are given by 

i( 2 ) 


Q2Q(m) 


r=3 |g| —r— 1 


8x a 


(P) 

"9192 Q „ 1 Z^ 

OX qi X q2 p=3 |h|=p 


EE 


^P0( m +2-p) 

dx h 


= 0, 


(65) 


( 66 ) 


^2^(771) ^+2 


W, 


( 2 ) 

Uj q2 


Qp~£(m+2—p) 


Z " =». (67) 

P=3 |ft,|=p 1 

where h and q are multi-indexes. The solution of equations (I63D - (1671) requires that the following 
normalization conditions are satisfied: 


-^J u< i m) (x)dx = o, ^ ^0W(x)dx = O, _i_^ T ( m )(x)dx = 0, 


( 68 ) 


where the £—periodic domain is the same defined in previous Section as C = [0, L] x [0, SL\. 

The averaged field equation (1551) . (1551) and m (or alternatively the sequence of PDEs (1551) - 
EZD), obtained by means of the proposed rigorous asymptotic procedure, are used in Sec. [|] of the 
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text of the paper for deriving the field equation of the first order (Cauchy) homogeneous continuum 
equivalent to the considered periodic thermodiffusive material. 

The approximation of the average field equations (I59l) - (l6l1l yielded by solution of homogenized 
differential problems of generic order m (1651) is more accurate with respect to that obtained by 
the assumption (|24l) . This implies also a more precise approximation of the solution of the micro¬ 
scopic field equation ([T )l -([8 Jl by means of the down-scaling relation (IT51) . (1TP1) and (1201) involving 
the m acros copic field (1521). As it is expla i ned fo r periodic elastic composites in lPeerlines and Fleckl 
( 2004 ) and iBacigalupo and Gambarotta ( 2012f) , the truncation of the average equations of infinite 
order (159D - (1G1I) at a generic order m with the aim to derive higher order field equations for gener¬ 
alized thermodiffusive continua may lead to problems in which the symmetries of the higher order 
elastic and thermodiffusive constants is not guaranteed. Moreover a loss of ellipticity of the govern¬ 
ing equation s can be observed. Asymptotic-variation al ho mogenization techniques simi l ar to those 
illustrated in Smvshlvaev and Cherednichenkol ( 2000l ) and iBacigalupo and Gambarottal ( 2012h rep¬ 
resent an appropriate and powerful tool in order to avoid these problems. The generalization of 
these methods to the case of elastic materials in presence of thermodiffusion is still missing in 
literature. 


B Symmetry and positive definiteness of elastic and thermod- 
iffusive tensors 


( 2 ) ( 2 ) ( 2 ) 

In this Appendix, the symmetry properties of the tensors of components n: iprn rj2 , rn qi , n , w qiq2 , and 
the ellipticity of the field equations (1001) and (1001) are demonstrated. 


B.l Symmetry and positive definiteness of tensor of components n^ qiq2 

(VS. Ciq 2pqi ) 

Let us consider the cell problem (121 H i . remembering that = 0, it becomes 

(<%ki N kL,i) ij +C% pqi j = 0, ( 69 ) 

where C^ kl are Q —periodic functions. The weak form of equation (1001) . using N^ q2 as Q —periodic 
test function, is given by 

((^<^ + ^ 0 ,; N lt) = °> ( 7 °) 

applying the divergence theorem to (1701) . and remembering that for the Q —periodicity of and 

2 the path integrals evaluated on the boundary of the unit cell Q vanish, we obtain: 

((q^SL+ c v P <n) N rlL) = (7i) 
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Using the result CD , expression (151[) -i can be written in the equivalent form: 


( 2 ) 

ipqi qi 


K( 


C E 4- C E 

w iq2pqi ' iq2kl iy kpqi,l 


) + ( 


,N, 


(i) 


c £ A- C e IX ' ■ 
^iq\pqi ' ^iqikl iy kpq 2 


-')) 


c £ A- C £ IX ~ - 
Ky iqipq 2 ' ^iqikl iy kpq 2 ,l 


. _L g 

rqii,J ' u rq 1 ^ij 


) 


1 (r? + c e /v (1) 4- /V m N (1) I pm A ax(i) , 

2 \ iq2pqi ' iq2kl kpqi,l ' y ijkl iy kpqi,l ' ijpqi J riq 2 ,j ' 

rym N w ■ sym \ N w ' 

-'rjkl ly kpq 2 ,l ' rjpq 2 J ly riqi,_ 

(c™ kl + 6 ri 6 jq2 + N^. + S rq2 S tJ ) ■ 

^ k P^ l 11 + ^kq lP ,l + ^ 

\ (c?jki 

(j^kpq 2 ,l + dkpfiq 2 l + ^kq 2P ,l + ^kq 2 ^lp^ 

as a consequence, we can observe that: 

jU 2 ) _ ifr'. i 

U ipqiq 2 ~ 2' ‘ lq2Pqi ' ‘ 

where the components Ci q2Pqi of the overall elastic tensor take the form: 

Ciq 2Pqi = — (c™ kl (j^riq 2 ,j + firidjq 2 + ,j + ^2%) ‘ (j^k Pqi ,l + ^kp^lq i + ^kq lP ,l + ^kqAlpj ^ i 

(74) 

Observing expression (1741) . it is easy to deduce that the tensor of components Ci q2Pqi is symmetric 
and positive definite. 


-iqipq2)i 


(72) 


(73) 


B.2 Symmetry and positive definiteness of tensors of components m q % 2 
and Wq% (vs. I< qiq2 and D qiq2 ) 

Remembering that tri;^ = 0, the cell problems E21) i . possesses the form 

+ 0 (75) 

where K™' are Q— periodic functions. The weak form of equation (1751) . using as Q— periodic 

test function, is given by 

( (KijM™ + K Z) . ) = 0 . ( 76 ) 

applying the divergence theorem to and remembering that for the Q— periodicity of If'" and 
Mql' 1 the path integrals evaluated on the boundary of the unit cell Q vanish, we obtain: 

( ( K Z<J + K Z) M S) = 0- (77) 
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Using the result (1771) . expression (15-11) mi can be written in the equivalent form: 


if 2 ) 


= t( (K m + K m M {1) ) + (K m + K m -M {1) ■) 

o \ v Q2Qi 1 qu qi,3 } ' 9192 1 fer ObV 


( 1 ) 

q2j^'"qi,j' 


K 


= -/K m + K m M (1) - + (K™M^■ + K m ^ M (1) ■ + 

2 \ 92 91 “ 9ir 92,J ' \ W 9i ,3 ' *9iy 92,* ' 

++*s) *&) 

+««)!+ 


9192 
1 ■ 
2 . 


= + s qiJ )) 


( 2 ) 

as a consequence, we can observe that rn qiq2 = K qiq2 , i.e. 


K, 


9192 


= + ^92)(< 1 i + <W) ■ 


(78) 


(79) 


Observing expression 02D, it is easy to deduce that the tensor of components K qiq2 is symmetric 
and positive definite. Since the the equations of heat and mass diffusion possess an identical form, 
the components of the tensors K qiq2 and D qiq2 have the same properties, and then the results 
obtained for the components of the overall heat conduction tensor can be extended to the case of 
the overall mass diffusion tensor of components D qi q2 . These components are given by the following 
expression: 


D 


9192 


( d ^(KI+ s ^Kj 



(80) 


C Overall elastic and thermodiffusive constants for bi-phase 
isotropic layered materials 


In this Appendix the explicit expressions for the overall elastic and thermodiffusive constant of a 
bi-phase layered material with isotropic phases are reported. The components of the overall elastic 
tensor take the form: 


Cim — 


-C 2 EgE b + C [{EgV b ) 2 - 2 EqVgEbVb + ( E b u a ) 2 - (Eg) 2 - {E b ) 2 ] - E a E b 

(C + 1)[C (E b ^g) 2 - E b ) + E a {v b ) 2 ~ Eg] 

(C + 1 )EgEb 


C 2222 — — 


C 1212 — 


Cl 122 — — 


c (E b {y a y - Eb) + E a {vb)* - Eg 
(C + 1 )EgE b 

2 [Eg + EgVb + C (EbVg + Eb)\ 

_ E a E b (v b + C "a) _ 

t(E b pa ) 2 - E b ) + E a {v b y - E a 


(81) 
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The components of the overall thermal dilatation tensor and diffusive expansion tensor are respec¬ 
tively given by 

Ana a — Bna a 
a n =- 7 -; 


0:22 = 


P 22 = 


C(E b (Pg ) 2 - E b )a a + E a ((v b ) 2 - l)a b 

(E b ((v a ) 2 - 1) + E a ((h) 2 ~ 1) 

0 A n P a ~BnP a 
pn --t-; 

C (E b (u a ) 2 - E b )p a + E a {Sy b ) 2 - l)/3 b 


(82) 


(E b ((ta) 2 - 1) + Ea((v b ) 2 - 1) 

where: 

A n = C 2 [E b (i>a ) 2 - E b ] + C [E a (v b ) 2 - E b v b + E b v b (v a ) 2 + E a v a - E a v a (v b ) 2 - E a ]\ 

Bn = C[ E a t> a (h ) 2 - E b + E b {v a ) 2 + E b u b - E b v b {v a ) 2 - E a v a ] + E a (h ) 2 - E a - 
All = (C + 1 )[C E b {("a ) 2 - 1 ) + E a ({h ) 2 - 1 )]. 

Finally, the components of the overall heat conduction and mass diffusion tensors become 


(83) 


An = 


D u = 


K b - (K a 

C + 1 ’ 

D b - (D a 


(C + l)K a K b 
K a +(K b ’ 

(C + 1 )D a D b 


K 22 = 


D 22 = 


(84) 

(85) 

( 86 ) 


C + 1 ’ “ D a +(D b ' 

D Down-scaling relations vs cells problems 

In this Appendix, we provide more details regarding the structure of the down-scaling relations 


Bensoussan et al. (1978); 

Bakhvalov and Panasenko 

(1984; Affairs ( 

1992); 

1993h; Meauid and Kalamkarov (1994) and 

Boutin 

(1996(1. the microscopic 


sented through an asymptotic expansion in the general form: 

+ OO 


Uh (*> 7) = J 2 £ ‘ u h ] (*> 7) = M i 0) (*- 7) + £u h ] (*> 7) + £2 4 2) (*> 7) + 
1=1 

+OO 

e (x, *)=E £ ‘° (l) (*> 7) = 0(o) (*> 7) + e0(1) (*> 7 ) + £20(2) (*, j ) + 


Z =1 

+00 


77 (x, y) = y; 7) = 7 ? (0) 7] + £?7 (1) j + e 2 ?/ 2) (x, * j + 


(87) 

( 88 ) 

(89) 


1=1 
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Substituting expressions (1571) . (1551) and (1551) into the microscopic field equations j7J), ([5]l. and re¬ 
membering the property £-f(x, £ = f) = (jj + §J£) ={§£; + ¥). . > we obtain 


£=x/e 


(C e 


+ e~ 


Cf 


ijhl 




dec; 


+c 


C| 


du 


(i) 


ijhl 


dxi 


+ u 


( 2 ) 

h,l 


dxj 


C e ( ( 1 ) 


+ . - (o,/") - (/3 |j?7 (0) ) 




c 


^ ( 9 u [ 2) (3)' 

yw i^r + ^ 


+ 


<9x 7 


J j 


(“‘A 2 ’),, - ST " (' 3 « ,(2, ) J - (' 3 «’ ?<1 ’) 


-2 


+£ 


+£ 


,J V dx, +t 

(dd^ 


k: 


i< 


( 2 ) 


d 

dxi 


ijhl \ dxi + h ’ 1 


~\~ bi — 0, i — 1,2, 


(90) 


? (!) 


9 r K e 

*r,j 


— f 

<9:r,- V 


I< 


dew 


lj ( dxj 


n(l) 


Kf 


lJ V a** 


+ 


3(3) 


d 

dx, 




■ + r = 0, 


(91) 


+ £ 


+£ 


+£ 


W), ( 


D e f^ + v W 

lJ \ dXj 


+ 


■L ( d « , 7) 


d 

l~ dxi 

d 

dxi 


D 6 (<>fL +r p 
10 \ dxj 

D-J&l + tf 


ZJ V dx, 


+ ■ 


■s = 0. 


(92) 


At the order e 2 from equation (l90l) we derive: 

{ C ijhl u hJ) . = ( 93 ) 

whereas from heat conduction and mass diffusion equations (1511) and (1551) we get respectively: 

(Kf/f) . = ff <°> (*), (Dtjvf) . = h {0) (x). (94) 

The interface conditions ©-CD, expressed with respect to £, become: 


,(0)t 


£esi 


= 0, 


C e 7 S°^n ■ 

°ijhl u h,l 


|esi 


= 0, 


(95) 
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il° {0) }] 

[b (0) ]] 


^esi 

|esi 


= 0, 


= 0, 




D v r f-j n > 


|gEi 

£gEi 



(96) 

(97) 


where Si is the representation of the interface E bewteen two different phases of the material in 
the non-dimensional space of the variable 
At the order e _1 equation (l90ll yields 


C « w l d Xl +Uh ’ 1 


+ ^ (Cf jhl u ( h °l) - (af/°>) - (^ i?7 (0) ) . = f?Hx), (98) 


at the same order, from equations and (j92]| we obtain: 

d 


' 00 (°) 

V ^7 






+ si( /r « e ,™) =9<i,<a;) - 

+ 7 K”. 1 ? 1 ) = ,, ' 1) w - 


and the interface conditions are given by 


,(i) T 


|esi 


= 0 ; 


( 0 ) 


re 9u h 

C ^ hl ~d^i 


+ u h,i - a ij e{0) ~ PijV 


(o) 


= 0, 


£eEi 


[[ 0 (1) ]] 


[[? 7 (1) ]] 


t =0; 

£g£i 

<H] 

t =0; 

K(^ )+ 9x,) n ‘I 


£e£i 


|esi 


= 0, 


= 0. 


At the order e°, the cells problems associate to equation ECT) assume the form: 

d 


C e . 

^ o' n 


du[ 1] (2) N 


ijhl \ dxi h ’ l 


dxj 


c: 


ijhl 


¥-2 


a 

dxj 


(a^6» (0) ) - (p ? ?? (1) ) - (/3^ (0) ) = /} 2) (*)j (104) 


whereas the cells problems correspoding to equations Eli) and (1921) are: 

d 


'86™ 

' i3 [~8^~ ' 


k;, (+ e [2) 


+ 


Df 


/ drp) 

V dxj 


( 2 ) 

■ 


dxi 

d 

dxi 


K*. 

lJ \ d Xj 


£>?■ W 1} 

lJ \ dxj + 


= 0 (2) (»), 
= h^ 2 \x), 


(99) 

( 100 ) 

( 101 ) 

( 102 ) 

(103) 


(105) 

(106) 
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and the interface conditions assume the form: 




£e£i 


= 0, 


CZjhl 


du 


(i) 


+ „g> - <,«(» - 


= 0, 




P (2) ]] 

[[?? (2) ]] 




|eEi 


= 0, 


= 0, 


D? 


? + 

aeK 1 ) ^ 

d Xj j 

1 rii 

|eEi 

= 0, 

(2) + 

drp^ N 

1 


= 0. 

,3 

Vxi j 

J 

£eEi 



(107) 

(108) 
(109) 


At the order e 2 , the solvibility conditions in the class of the functions Q —periodic with respect 
to the fast variable $ implies that f- 0> (x) = g(°\x) = h^°\x) = 0, then the cell problems (1M1) - (I91I) 
become: 

. = o, (Kt 3 ef ). = o, (pijrif). = o, (no) 

as a consequence, the solution of problems (II101) does not depend by the fast variable £ and then 
«l 0) (a:,O = U h (x), 0 (o) (x, £) = Q(x) and r/ 0) (a;,£) = T(x). 

At the order e _1 , the solvability conditions in the class of the functions Q —periodic with respect 
to the fast variable £ together with the interface conditions (I101I) - (1103D yield to 


=/" w. 

= S m M, = h m {x). 


dxj 


(HI) 

( 112 ) 


For the Q— periodicity of the functions C-j hl , a ?•, /3fj, Kf- and Df-, we have f-- l> (x) = g^(x) = 
h W (a:) = 0, and then at this order the solution of the fields equations assumes the form: 


(€)©(*) + N^m(x) 


— a d 1 ) 


dU , 


(i)/ 


r(!)/ 


(113) 


0 w (x,£) = m<U(s) 


<90 

dxri 


v w M = w£>{t) 


or 

’ dx„ 


, (H4) 

where , N^\N^\ and are the same fluctuations functions introduced in Section[3] 
Substituting expressions (I113I) - (I114I) into the cell problems (I111I) - (I112I) and considering the interface 
conditions (110111 - 1)10311 . we derive: 




C £ = 0 
ijpqij ’ 




— Q/?. . = 0 
LX ZJJ W, 




+K s =0 fD E W (1)N ) +1F £ =0 

*3 ?lj^ T m,* ’ *3 9 l Jy i T rv iqi,i u ’ 


and then the interface conditions (I101I) - (I103I) become; 


[[IV; 


(i) ] 

hpq if 


£eSi 


= 0 , 


QjTiZ (j^hpqi.,1 4" fihpfilq^J Tlj 




= 0 , 


3fi,i = 0, (115) 
(116) 

(117) 
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IWA 


a: 


IN, 

K’l 

[[< 1} ]] 


|e£i 


|eEi 




|esi 


= 0, 


= 0 , 


= 0, 


= 0 , 


Kfj [Mq2i ^91 i) ni 

Dfj (Wqjj + 5 qi j) rii 


|esi 

^€Ei 

£g£i 


= 0 , 
= 0 , 
= 0 , 
= 0 . 


(118) 

(119) 

( 120 ) 
( 121 ) 


The solution of the cell problems (II 151) - (II16D taking into account the interface conditions (II 1 71) - (112 ID 
provides the Q —periodic perturbation functions N^J qi , N^\N^\ and W q ] 1 . 

Taking into account the solvability conditions in the class of the functions Q —periodic with 
respect to the fast variable £ and the interface conditions (I107I) - (I109I) . the cell problems (I104D - (I106D 
associate to the e° become: 


(Cf qipl 

+ CilhjNhJ qit j ) 

d 2 U p 

(QmAj 


M dx t 

= /f >w, 

( 122 ) 

dx qi dxi 

{ K » 

("£+■*»)} 

d 2 Q 

= g (2) {x ), 

K(<j+«»)) 

d 2 T 

= h^ 2 \x). 

(123) 

dxidx qi 

dxidx qi 


these cell problems possess a solution satisfying the conditions (I107D - (I109D in the form: 

,( 2 ) („ e\ - at ( 2 ) d 2 £/ p ^ a7(2) _9B_ A y(2) dT 

Vqi 


i ( > (x£)=N'- ) p _I N ( 1 — I N ( ’ 

h 1 >41 h pqiqi g Xqi d Xao + h ^dx„A hqi dx, 

d 2 e 


0W(x,Z) = m£1 


qiq2 dx q2 dx qi 


rj( 2 \x 


£) = tr (2) 

sy vv qiq 


Q l 

d 2 T 


qiq2 dx q2 dx qi 


(124) 

(125) 


(2) ~ (2) A (2) (2) (2) 

where N h ' q2 , N hqi N hqi , M qiq2 and W qi ; l2 are second order fluctuations functions just introduced 
in Section |3 As a consequence, from the cell problems (I122I) - (I123D we derive: 


C e 

^tjhl ly hpq 1 q 2 ,l 


C e N^ 

^ijhl^hqx 


)j + ( 


C e ^ 

tjhq 2 hpq x J 


■CL^+Ct q2hj NlJ qi J = (C; 


- y iqipq2 


c, 


ijhq 1 I 'h 


N 


(i) 


). 


+ ^iqikjN^j - ( CX 


0 


\ N K } 

iqipq2 ' ^hj 1 ' hpqx, 


- < = {Ct qihj m - a 


. + ( C !j hqi N™) _ + Ct qikj Nfj - (^.<)) _. - Pi qi = (c? qihj N% - (3. 

. + ( A % 2 A/ 4 1} ) . + A %. (Kj + <W) = (Kj (<J + <W) ) > 


(d s W < ' 2) + (D e 1T (1)N ) + D e (w^-+6i„\ = /D e ■ ( 

\ v 9192 ,j) i “ \*92 vv qi J iq 2 \ vr qi,j ~ u JQi J \ q 2 j \ 
and then the interface conditions (I107I) - (I109I) become; 

(^•«< ) gi 924+^92<y% 


lwW+6 jqi 


q i ,j 


))• 




(2) l 

hpqiq 2 \ 


|eSi 


= 0, 


£e£i 


= 0, 


), (126) 

)> (127) 
}, (128) 
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(130) 
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[[iv/ 2) 


hqi 

[[^ 


£e£i 


= 0 , 


[ C W*&,i + C Zjh q M ) -< 


£esi 


= 0 , 


( 2 ) 

hqi 


|eSi 

mZt 

[[w'SJ 


= 0 , 




£esi 


= 0 , 


= 0 , 


Kf, M. 


( 2 ) 




?192 JJ ^ gEi 

The solution of the cell problems (11261) - 


S'f, IT, 


( 2 ) 

9192 ,j 


^GSi 

|esi 


|eEi 
= 0 , 

= 0 . 


= 0 , 


(132) 

(133) 

(134) 

(135) 


K ] ) n > 

KP-fijWW) 

+ M qi^ ^JQ2 ) TH 

W to S 3<to) fli 

taking into account the interface conditions (I131I) - (I135I) 
provides the Q —periodic perturbation functions , N^, N^, Mq 2 \ 2 and Wq 2 \ 2 . 

The general procedure here reported can be applied to higher order cell problems for deriving 
the averaged field equations of infinite order, which assume the form: 

(C iq lP l + C'il.h.-iNj' ^ 

r(i) 


V &* U P , C e jjW _ n eK^__, C e ftW _ ae\&L 

iq!pl T '^zlhj 1 ' l hpqi,j/Q x ^ Q Xi \ ( ~‘llhj I ' l h,j ^"dxi ^ ll ' dXl 

d 3 U v 


'ilhj ly h,j 

HD 


+(C,U,<>, + C t,,uKL,„l) g^tsZ, + ~ <Z M k ) )g^ 

d 2 T 


+ {CL 1 ha 2 fth + CLhlN&A - 0LW&) 


( K Li ( 


Mqi'j + Sjq 1 


iq2hl iy hqi,l 

(i) 

7l,. 

.+ r(a:) = 0 

( D ( 


iq2 91 dx qi dx q2 


+ bi(x) = 0, 


<9 2 e 

qi' JJj q 2 

(136) 


—-h (K e M (1) + I< £ M (2) \-— 

dx q2 dx qi \ 9392 92 9 « 9i92 ’ 9 / dx n ,dx no dx. 




<9 2 T 

dxq 2 dx qi 


+ (D e W (1) + D E W 

' \ 9392 92 ' q3j c 


v qi v<vq 2 i. 


d 3 T 


93 


( 2 ) 

9192.J/ dx qi dXq 2 dx q3 


+ s(a:) = 0. 


(137) 


(138) 


Note that applying the permutation of the saturated indexes, (I136D . (11371) and (11381) become iden¬ 
tical to the averaged field equations derived in Appendix [A] The structure of the down-scaling 
relations is defined by the solutions of the various cells problems associate to the different orders of 
the asymptotic expansion. These down-scaling relations assume the form: 


Uh 


(«■=) 


[4(*) + e (AC(€) 


dU p (x) 


dx +Nl 1 \t)Q(x)+N^((i) T(*)) + 

, -2 « * r ( 2 ) ( ^d 2 U P {x) A ~ T ( 2 ),„dO(x) A y( 2 ). ^ 0T(*) , 

+ £ (^P9i92«)^-^-+^ l (0^— + N h qi (0^— ) + 


tWi 


J^=CC/ E 


< x \ 


hi) 

— 


0(x)+ £ mW(O 


90(x) , 2 %r( 2 ) d 2 Q(x) 


Sin 




dXq-^ 0Xq 2 




- £=Z/e 


(139) 

(140) 

(141) 


The relations (11391) . (1140p and (11411) are identical to expressions (fl2l) . (fl3l) and (IT4)) introduced in 
Section [3] and used for developing the homogenization method illustrated in the paper. 


40 



































































